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Abstract 

We study the T odd correlations induced by CP violating anomalous top-quark couplings at 
both production and decay level in the process gg — > tt — > (bfi + u^)(bfj,~ v^) . We consider several 
counting asymmetries at the parton level and find the ones with the most sensitivity to each of 
these anomalous couplings at the LHC. 
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I. INTRODUCTION 



With the upcoming start of the LHC in mind, we consider the possibility of looking for 
CP violation in high energy processes involving the production and decay of top quarks. 
There are many observables suitable for this purpose that have been studied before in the 
literature 

Qaaaaaaai. 

In this paper, we consider T odd triple product correlations 
of the sort first discussed in Ref . jl] . 

We keep our study as model independent as possible by considering a scenario in which 
only the standard model (SM) particles are relevant and the CP violation is parametrized 
by anomalous top-quark couplings affecting both the production and decay vertices 0] • The 
ti production process is modified relative to the SM by the interaction 

Ccdm = -ig s ^icr^ v 'j 5 G tlu t, (1) 

where g s is the strong coupling constant and G^ v is the usual gluon field strength tensor. This 
interaction results in a modification to the tig vertex as well as in an additional, "seagull" , 
tigg vertex. For the decay process t — > bW + we write the most general decay vertex (and 
the corresponding one for i decay), 

T " Wtb = ~7I Vtl U{Pb) ^ {flLPL + f * PR) ~ ^ {Pt " P^ Pl + ^)] (2) 



This vertex can be derived from a dimension five effective Lagrangian as in Ref. [10j, but 
unlike the case of Eq. [TJ the effective Lagrangian does not generate other vertices that affect 
this calculation. For the remainder of this paper we will use Vtb = 1, f± = 1, fx = 
and fi = as in the SM, and allow for new physics only through the coupling f% which 
is the only one that can interfere with the SM to produce T-odd correlations. To generate 
T-odd observables the coupling _/vf must have a phase but this phase does not have to be 
CP violating. We thus write f% = /expi(0y + Sf) using (fif to parametrize a CP violating 
phase due to new physics and 5f a CP conserving phase arising from real intermediate states 
at the loop level. Some of these CP violating couplings have also been studied in the context 



of a future linear collider 11 



Our study in this paper corresponds to the numerical implementation of the results 
presented in Ref. 12], where comparisons to the previous literature were made 0, 0, S]- 



We start from Eqs. 10 and 32 of Ref. [12j which give the spin and color averaged matrix 
element squared containing the T odd triple product correlations for the process gg —>■ it — > 



b/j, + i>nb{i fu. For convenience we collect the results of Ref. 12| in the Appendix. Numerical 



studies similar to ours have been carried out before for a subset of the observables we consider 



here, and we compare our results to the most recent ones in the literature 13j, ll4j. The 



numerical studies are performed with the aid of MADGRAPH [15|, [16|, [1 



II. OBSERVABLES 



In Ref. 12j three independent correlations were identified to test CP violation in tt 
production. Here we rewrite them with a small modification: we have reshuffled factors of 
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t — u between the form factors and the correlations to ensure that both the form factor and 
the correlation are even under the interchange of the two initial protons. We concentrate 
initially on W decaying into muons, for which the T-odd correlations are 1 : 

Ox = e(p t ,pt,p^ l +,p tl -) 

2 = (t-u)e(p lt +,p li -,P,q) 

£>3 = (t-u) (P ■p^+e(p lx -,p t ,pi,q) + P •Pp-e(Pp+,Pt,Pf,q)) (3) 

with q = pi — P2 and P = pi + P2 being the difference and sum of the incoming parton mo- 
menta. The spin and color averaged matrix element squared that contains these correlations 
is given by, 

\M\ 2 CP = C 1 (s,t,u)0 1 + C 2 (s,t,u)0 2 + C 3 (s,t,u)0 3 , (4) 

where the form factors Ci.2,3 were computed in Ref. and we reproduce them in the 
appendix for convenience. 

We begin by studying each of the three terms in Eq. H] separately, considering the lab 
frame distributions da/dOi for the three correlations. In each case we isolate the CP odd 
form factor C, by constructing the integrated counting asymmetry 

a _ N events (Oi > 0) — N events (Oi < 0) 

' = N events (O t > 0) + N events (Oi < 0) ' 1 j 

The observables used to construct the Ai are not realistic in that not all the momenta 
appearing in them can be reconstructed. To address this issue we replace those observables 
assuming that for each event it is only possible to reconstruct the momenta of the two muons 
/r^, the two b,b jets, and the beam direction. The correlations under this assumption can 
be obtained from Eq. [3] with the substitutions 

Pt^ Pb+ Pu+ Pi ->Pb+ Pa- 

P -> Pb+Pn+ +Pb+Pfi- q^q = Pi- P 2 - (6) 

We have defined a four-vector q, as the difference between the two beam four-momenta. The 
factor t — u could get modified by writing it as (t — u) = q • {pi — Pt) with the substitutions 
implied by Eq. [61 However, all one needs is a factor linear in q so we choose the simpler 
form (t-u) q - (p n - ~Pn+)- 

All this results in the correlations O, 

01 = t(Pb,PhPu+,Pu~) 

02 = q- (Pu+ -p»-)t(p»+,P»-,Pb+Pt,q) 

03 = q - (Pa+ - Pa- )e(Pb,Ph~, Pa+ + Pa- ,q), (7) 

and their associated counting asymmetries A*. It is easy to see that the correlation O3 gives 
rise to both O2 and O3. From the experimental perspective, 2 is most desirable as it is 
the only one that does not require distinguishing between the h and b jets. 



Here we use the Levi-Civita tensor contracted with four vectors e(a, b, c, d) = e^apa^b^c 01 ^ with the sign 
convention 60123 = 1- We also use s, t, u to refer to the parton level Mandelstam variables for gg — > it. 
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When CP violation occurs in the decay vertex, the spin and color averaged matrix element 
squared containing the T-odd correlations was written in Ref. 12j as 2 



\M\t = fsm(<f) f + 5f)e(p t ,p b ,p e +,Q t ) + fsin(<f) f -Sf)e(pt,pi,pe-,Qi)- (8) 

All the terms in Eq. [8] contain three four-momenta from one of the decay vertices so the 
correlations 01,2,3 defined previously may not be the best to measure these couplings. Guided 
by the form of Eq. [HI we define the following three correlations for this purpose 

04 = e(P,p b - Pi,Pn+,Pn-) 

05 = e(Pt,Pi,Pb+Pb,P»+ -Pn~) 

6 = (t-n)e(P,p 6 + ps,p M + -p^,q)- (9) 

These three correlations explicitly show that in order to test for CP violation in the decay 
vertex it is necessary to compare the decay of the top quark with that of the anti-top quark. 
This is accomplished in Eq. M with the use of the linear combinations of pb ± pi as well as 
p M + —Pfj,-. These constructions are CP-odd, and as such they isolate the phase 0/ in Eq. 
Later on we discuss alternative constructions to isolate the phase 5f. Interestingly, after we 
use the substitutions of Eq. [6] to account for the fact that the top four-momenta cannot be 
reconstructed completely, no new correlations are needed: both 04,5 become proportional 
to 0i and 6 becomes proportional to 2 . 



III. NUMERICAL ANALYSIS 



We start from the standard model process gg — > tt — > b^Vyba z/ M implemented in MAD 



GRAPH according to the decay chain feature described in Ref. 17| . This decay chain feature 



is chosen for consistency with the approximations in the analytical calculation of the CP 



violating interference term presented in Ref. 12J , in which the narrow width approximation 



is used for the intermediate top quark and W boson states. The expressions from Ref. [12j 
are then added to the spin and color averaged matrix element squared for the SM (which 
MADGRAPH calculates automatically) and the resulting code is used to generate events. 
This code is then missing the terms that are completely due to new physics: those propor- 
tional to the anomalous couplings squared. This approximation is justified because those 
terms do not generate T-odd correlations. In addition, as long as the conditions that allow 
us to write the new physics in terms of anomalous couplings remain valid, their contribution 
to the total cross-section is small. 

For event generation we use the default MADGRAPH cuts requiring the top quark and 
W boson intermediate states to be within 15 widths of their mass shell, the pt of both muons 
to be larger than 10 GeV and 77^ < 2.5. We also use SM parameter values as in Madgraph, 
except for m;, = 0; and we use the CTEQ-6L1 parton distribution functions. Furthermore, 
we only include the gluon fusion initiated parton processes. This procedure leads to a 
cross-section a{pp — > tt — > b^v^^T Vy) « 4.3 pb, which is independent of the anomalous 



Note that there is a typo in Ref. 12| where cj>f and 5f are reversed. 
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top-quark couplings d, f. The total cross-section does not depend on d, f because we are 
only keeping deviations from the SM that are T odd. These terms are also P odd, and as 
such they integrate to zero. 

The 4.3 pb cross-section we obtain underestimates the known theoretical cross-section 
by almost a factor of two. The latter is derived by multiplying the next-to leading order 
(NLO) total cross-section for tt production at the LHC with the W — > fiu branching ratios: 
833 pb x (0.1057) 2 « 9.3 pb. About 20% of this difference arises from the NLO corrections 
and most of the rest from our decay chain approximation. We expect the asymmetries 
we compute to change slightly when the decay chain approximation is relaxed and NLO 
corrections are included. Numerically we have checked that the asymmetry increases with 
the cross-section as we relax the decay chain approximation. We cannot quantify rigorously 
this increase because our starting (analytic) formula was computed in the narrow width 
approximation. In addition, since the asymmetry is an interference effect between the new 
physics and the SM amplitude, we expect it will get modified by about 10% by NLO QCD 
corrections. 

We first estimate the counting asymmetries of Eq. [5] by generating 10 6 events for each of 
the four cases: d = 5 x 10~ 4 GeV -1 ; /sin^ = 5 x 10" 4 GeV" 1 ; / sin 5/ = 5 x 10~ 4 GeV" 1 
and d = f = 0. These cases correspond to CP violation in the production vertex, CP 
violation in the decay vertex, strong phases in the decay vertex and the lowest order SM 
respectively. According to our discussion in Appendix [HI we require that the signals be 
larger than the 3a level for statistical fluctuations, ~ 3 x 10 -3 . Asymmetries found to 
be at this level or below, are estimated once again by increasing d or / by factors of ten, to 
5 x 10~ 3 GeV -1 . If the results scale linearly with d (or /), they are treated as true CP-odd 
asymmetries and entered into Table [H Otherwise they are treated as zero and examined 
further in Appendix [B] by increasing the number of generated events and/or the size of the 
coupling. All the zero entries in Table [I] satisfy our a-priori expectations. All the observables 
0i_6 are truly CP odd, so they do not get contributions from the SM (d, f = 0) nor from 
the final state interactions (sin0j = but sin 5/ ^ 0). The case of Ai_ 3 as induced by 
fsincpf is not clear cut as seen in Appendix [Bl but we do not pursue it further since there 
are larger, more promising, asymmetries to constrain this coupling. 
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TABLE I: Integrated asymmetries for d, fsincpf or /sin<5y = 5 x 10 4 GeV 1 and for the SM 
(f,d = 0). 

We now consider two types of effects that can dilute our estimated asymmetries. First, 
since the neutrinos are not detected, full reconstruction of the top four-momentum is not 
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possible. As discussed above we deal with this complication by considering the asymmetries 
Ai of Eq. [7] instead of the Ai. Second, we want to introduce cuts to suppress known back- 
ground and see their effect on the asymmetries. Specifically we know that the cross-section 
for pp — > bbfi + fi~X (no intermediate tt) within the SM is significantly larger. This process 
is dominated by strong production of bb pairs with the lepton pair attached to a photon or 
Z boson. With the same default Madgraph parameters and cuts, we estimate this cross- 
section to be about 24 pb, or some 5.6 times larger than that for the process of interest to 
us. This SM process is CP conserving so it cannot fake any of our CP odd asymmetries. 
Nevertheless, it will make any non-zero Ai smaller because it increases the total number of 
events. To reduce this background, we introduce two standard sets of cuts. The first is a 
set of acceptance and separation cuts for the muons and b quarks are: 

Pt{^) > 20 GeV pr(b,b) > 25 GeV 

\r](b,b, ^)\ < 2.5 AR(bb) > 0.4. (10) 

This set of cuts already reduces the background cross-section by a factor of 20, to about 
1.2 pb. A second set of cuts includes Eq. [TOl and also requires missing transverse energy, 
P T , to select events originating from semi-leptonic t, i decays, 

$ T > 30 GeV. (11) 

The requirement of a minimum $ T completely eliminates the background in our parton level 
study, although this will no longer be the case when detector effects, initial state radiation 
and final state radiation are included. 

In Table HT1 we show the effect of these cuts in our asymmetries. Comparing Tables HI and HT1 
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TABLE II: Integrated asymmetries for d or fs'mcpf = 5 x 10 4 GeV 1 with the cuts defined in 

Eqs. ma nil 

we see only small changes, with some asymmetries actually being larger after the cuts are 
imposed. This, of course, is due to a smaller denominator in the respective Ai reflecting 
the smaller number of events. Within our approximations, the cross-section becomes 2.6 pb 
with the cuts of Eq. QJJ an d 2.3 pb with both sets of cuts indicating a loss of 40 — 50% in 
the total number of events. 

Next we calculate the asymmetries associated with the correlations Oi^ (which do not 
depend on full reconstruction of the top four- momentum) and show these results in Table [TTT1 
For d — 5 x 10~ 4 GeV -1 and /sin^j = 5 x 10~ 4 GeV^ 1 we present results corresponding 
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TABLE III: Integrated asymmetries without full top momentum reconstruction for d or / sin <j)f 
= 5 x 10" 4 GeY" 1 with the cuts defined in Eqs. HQ| ED 
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FIG. 1: da/dOi and da/d&x distributions for the cases d = (SM) and d = 5 x 10" 4 GeV -1 as 
well as d<j/d& 2 for / sin (f> f = 5 x 1(T 4 GeV" 1 . 

to the two sets of cuts. The largest asymmetry for a d coupling is A\ which requires 
distinguishing between the b and b jets. Interestingly the most promising asymmetry for 
probing CP violation in the decay vertex is A 2 , which does not require distinguishing the b 
and b jets. To gain more insight into these asymmetries we show the differential distributions 
da/dOi and da/d&x for d = 5 x 1(T 4 GeV" 1 as well as da/d0 2 for / sin cj) f = 5 x 10" 4 GeV" 1 
in Figure [H 

The distributions shown in Figure [1] can also be used to estimate the sensitivity of LHC 
to the anomalous couplings as was done for the case of d and 0\ in Ref . [l3| . In this paper 
we rely exclusively on the counting asymmetries since we are ignoring detector and other 
issues that are more important for the distributions. We summarize our results for the most 
promising asymmetries in terms of the dimensionless anomalous couplings 



d t = dm t , f t = fm t 



(12) 
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for m t = 171.2 GeV as follows. The largest asymmetry is 

A x = 1.17 d t . (13) 

Measurement of this asymmetry requires reconstruction of the top pair four-momenta. As- 
suming that only muon and b four-momenta can be reconstructed, our asymmetries are 

A x = 0.64 d t - 0.072 f t sin (j) f 
A 2 = -0.041 dt - 0.32 / t sin 0/ 

A 3 = 0.21^ + 0.047/^^0/. (14) 

We now estimate the 5a sensitivity of LHC with an integrated luminosity of 10 fb^ 1 to 
the anomalous couplings. Using Eq. IB1I with iV = 23A; events per year (corresponding to 
a = 2.3 pb) we see that 5a sensitivity requires A{ > 0.033. From A\ and A<i respectively, 
we find setting only one anomalous coupling to be non-zero at a time, 

\dt\ > 0.05, \d\ > 3.0 x 10" 4 GeV" 1 (15) 
\f t sm(p f \ > 0.10, \fsm<p f \ > 6.0 x 10" 4 GeV" 1 

The study of Ref. [ill] found that the LHC with an integrated luminosity of 10 fb^ 1 could 
achieve a 5a sensitivity to d t > 0.05 by considering an observable proportional to A\. The 



agreement with our estimate is very good although the analysis in Ref. 13J is much more 



complete. The conclusion by Sjolin in Ref. 13j is reached by analyzing both the counting 



asymmetry and the distribution proportional to da jdO\ including both di-lepton and lepton 
plus jets channels from the W + and the W~ in the semileptonic t and i decays. His study 
also includes many background and reconstruction issues that we have completely ignored. 
However, we can roughly compare the two results: we have used only the channel with fi + 
and /i~ from the W + and W~ decays and have looked only at the counting asymmetry. The 



comparable 5a sensitivity result of Ref. [13j , from the counting asymmetry in the di-lepton 
channel is based on 13k events that survive all the cuts. After our less restrictive cuts, our 
5a sensitivity result is based on 23fc events resulting in about 30% more reach. We conclude 
that the two estimates are consistent with each other. 



IV. STRONG INTERACTION PHASES 

In the previous sections we have constructed observables designed to isolate true CP 
violation from the T -odd triple products. However in some cases it may be desirable to 
isolate CP conserving but T-odd effects induced by strong interaction (unitarity) phases 
and we address this issue now. We begin by choosing two observables to construct counting 
asymmetries. From Table [TTT1 it is clear that O2 is the most sensitive observable to phases in 
the decay vertex. We can easily construct a CP-even version of this observable by replacing 
the sum of b and b jet momenta with their difference. Since it is hard to distinguish these 
two jets, we would like to have at least one example of a CP even observable that does not 
require this separation. Thus we consider 

O a = q- (?V +P»-)e(p»+,P f M-,Pb + PhQ) 

o b = q-(p»+ -p»-)e(p»+,P»-,Pb-Pb,q)- ( 16 ) 



S 



Proceeding as in the previous section we find results shown in Table IIVI 
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Eq.[TD] 
Eqs. nm HI 



TABLE IV: Integrated CP-even asymmetries with / sin Sf = 5 x 10 4 GeV 1 with the cuts defined 
in Eqs. QZ2 dU 

Using the larger of these two asymmetries we can write 

A b = -0.32^ sin 5 /; (17) 

from which we conclude that the LHC with 10 fb _1 will have a 5a sensitivity to 

\f t sin 5 f \ > 0.10 |/ sin 5/1 > 6.0 x KT 4 GeV" 1 . (18) 

The magnetic transition factor / is induced by one-loop QCD corrections at the level / = 
8 x 10~ 5 GeV -1 Lsl ]. and it has been claimed that this level can be reached at LHC 10]. The 
la contraint for 10 fb" 1 has been estimated at -2.6 x 10" 4 GeV -1 < / < 1.4 x 10" 4 GeV^ 1 
by setting other anomalous couplings to zero in Ref. 14|. This constraint was obtained using 
T-even observables and is somewhat better than what we find with the T-odd correlation 
O a . Of course, our T-odd observables require a non-zero absorptive phase 5f, unlike the 
observables used in the study of Ref. 



14 



Although it is easy to see that QCD corrections 
to the tree-level weak vertex can introduce such a phase, it has not been calculated yet to 
our knowledge. 



V. SUMMARY AND CONCLUSION 



We have presented the results of a numerical implementation of the results in Ref. |12 
using MADGRAPH for event generation at the parton level. We have used these events to 
estimate the sensitivity of the LHC to CP violating anomalous top-quark couplings. 

For the case of the coupling d t that parametrizes CP violation in the ti production 
process, we find that the LHC with 10 fb _1 can rule out values larger than 0.05 at the 5a 
level. We reach this conclusion by considering only counting asymmetries in the di-muon 
channel. Specifically, our result follows from A\ which is the most sensitive asymmetry 
that does not require complete reconstruction of the t and i four-momenta. Our result is 



consistent with the estimates of Sjolin in Ref. [13J. Further refinements to our study that 
can be carried out by the experimental collaborations along the lines of Ref. 13] include: 
going beyond the parton level to include detector effects; using lepton+jets channels; using 
distribution shapes in addition to the counting asymmetries. Our paper has shown that there 



are other asymmetries, beyond the one studied in Ref. 13], that can provide comparable 



levels of sensitivity, and/or different handles on the analysis. For example the correlation 
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Ai used by Sjolin requires distinguishing between the b and b jets. Although this may be 



possible, as discussed in Ref. 13j, we have presented an alternative observable, A 2 , where 
this is not necessary. 

For the case of the anomalous coupling f t sm<f)f that parametrizes CP violation in the 
decay vertex, we find that the LHC with 10 fb^ 1 can rule out values larger than 0.10 at the 
5cr level. This result is based on the counting asymmetry A~i- Using Ref. [13| as a guide, we 
may expect a loss of a factor 2-3 in sensitivity when going from our parton level study to a 
detector level one. 

For the case of CP even but T odd correlations induced by strong phases in top decay, 
we find that the LHC with 10 fb _1 can test values as small as f t sin 5f = 0.10 at the 5a level. 
This estimate is consistent with that of Ref. 14] for f t using different observables. 

We have shown that, at least at the parton level, it is possible to significantly reduce 
or even eliminate the dominant background with a simple set of cuts that have a relatively 
small effect on the signal. CP conserving background can arise from the fact that the initial 
state at LHC (pp) is not a. CP eigenstate. This issue has not been addressed in detail so far, 
but a simple estimate in Ref. [6j puts it several orders of magnitude below the sensitivity 
that we have found. 
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APPENDIX A: FORM FACTORS 



The form factors for CP violation in the production vertex appear in Ref. [12j split 
into three different contributions. The three contributions correspond to those of different 
gg — > ti diagrams and were split that way to facilitate checking them. In this paper we are 
only interested in the total result for use in numerical routines. This corresponds to Eq. 16 
of Ref. [12| which we write here explicitly: 

C X = 6s2 J^l uW (9(* - uf - 2s\t - uf - 7s\t - uf (Al) 

+32m 4 (7s 4 + 9s 2 (t - uf) + 4m t 2 (7s 5 + 7(t - ufs 3 + 18(t - ufs)) 

C * = ^gffi (_ 7 s 4 - 2(t - ufs 2 + 9(f - uf + Am 2 (23s 3 + 9(t - ufs)) 

bs z (s z — (t — uyy 



and C2 = §£3. Notice that these form factors differ form those defined in Ref. [12J by 
factors of t — u according to Eq. [31 

The Qt and Qi that appear in Eq. [8] are linear combinations of available momenta that 



act as spin analyzers for the t and t respectively. They were also given in Ref. 12j separately 



for three different contributions from the different diagrams in gg — > tt. Once again, we are 
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only interested in the sum of these contributions in this paper so we present it here explicitly 
for convenience: 



Qt = K U fe2(s2 _ ( ;_ M)2)2 { ("32(7s 4 + 9(t - M ) 2 s 2 R + 8s(s 2 - (t - uf) 

■(Is 2 + 27(t - uf)m 2 + (29s 2 - 45(t - uffs 2 - (t - u) 2 ) 2 ) p e - 

+2 (-s(29s 4 - 74(t - u) 2 s 2 + 45(t - w) 4 + 4m 2 (7s 3 + 9(t - u) 2 s))^- • p t 

+s{9{t - uf + 2s(54m t 2 + 17s) (t - uf + (84m 2 - A3s)s 3 )p e - ■ p- t 

+ (t- u)(A5(t - uf + 2(54m 2 - 37s)s(t - u) 2 + s 3 (29s - 44m t 2 ))^- • q) p- t 

+2 ((t - u)(-45(t - u) 4 - 2(54m t 2 - 37s)s(t - u) 2 + (44m 2 - 29s)s 3 )p t - ■ p t 

+(t - u)(-45(t - m) 4 - 2(54m 2 - 55s)s(t - u) 2 + (44m 2 - 65s)s 3 )^- • p f 

+ s(45(t - w) 4 + 2(18m t 2 - 37s)s(t - w) 2 + s 3 (28m 2 + 29s))p € - (A2) 



Qi = 6g2(g2 _";_ M)2)2 ( (-32(7s 4 + 9(f - u) 2 s 2 )mt + 8s(s 2 - (t - u) 2 ) 

•(7s 2 + 27(* - uf)m 2 + (29s 2 - 45(t - w) 2 )(s 2 - (t - w) 2 ) 2 ) p e+ 

+2 (-s(29s 4 - 74(t - u) 2 s 2 + 45(t - u) 4 + 4m 2 (7s 3 + 9(t - w) 2 s))p £+ ■ p f 

+s(9(t - w) 4 + 2s(54m 2 + 17s) (t - u) 2 + (84m 2 - 43s)s 3 )^+ • p t 

- (t - u)(45(t - uf + 2(54m 2 - 37s)s(t - u) 2 + s 3 (29s - AAm 2 t ))p l+ ■ q) p t 

+2 (-(t - u)(-45(t - u) 4 - 2(54m 2 - 37s)s(t - u) 2 + (44m t 2 - 29s)s 3 )^+ ■ p { 

-(t - u)(-45(t - uf - 2(54m 2 - 55s)s(t - uf + (44m 2 - 65s)s 3 )^+ • p t 

+ s(45(t - uf + 2(18m f 2 - 37s)s(t - uf + s 3 (28m 2 + 29s))p £+ • q) q} . (A3) 



The factor Ku defined in Ref. [12f] is given by: 



K u = 16 (n 2 a 2 s g 8 ) (p b ■ p u ) (p~ b ■ p D ) 



7T \ / 7T 



m t T t J \M W Y W 

x 5(p 2 - m 2 )5(p 2 t ~ m 2 t )5(Pw + ~ M 2 w )5{p 2 w - - M 2 W ). (A4) 

For our numerical implementation we rewrite all four delta functions as the respective Breit- 
Wigner distributions behind them, for example: 



7T 

mttj"^ (p 2 t -m 2 f + T 2 m 2 



S(p 2 - m 2 t ) -> — 2 (A5) 



APPENDIX B: SMALL ASYMMETRIES AND STATISTICAL FLUCTUATIONS 

For our numerical analysis of the counting asymmetries in Eq. E] it is important to know 
the level at which these can be faked by statistical fluctuations. For a CP conserving 
background (such as the SM) which does not generate the asymmetry, the probability for a 
given event to fall in the bin with Oi > is p = 1/2. The number of events falling in this 
bin out of a total number of events N, follows a binomial distribution with mean Np and 
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standard deviation a/ Np(l — p). It follows that N events due to a CP conserving process 
lead to the expectation 

Ai = 0±(-^=] . (Bl) 

Several asymmetries were found to be zero in this paper and we now show that those 
asymmetries are really zero and not just numerically small. To do this we show results 
obtained by varying the number of generated events and/or the size of the anomalous cou- 
plings to demonstrate that certain small asymmetries are consistent with expectations from 
statistical fluctuations, and not with CP violation. For the SM case we show in Table IVl the 
asymmetries that we obtain when generating 10 6 and 10 7 events with MADGRAPH using 
the cuts in Eqs. [TUJ [TTJ For comparison we show the 3<r statistical error from Eq. IB II above. 
The results in Table [V] show that at the 99% confidence level the SM does not induce any 



N 


3 


A x 


A 2 


^3 


A± 


A 5 


A 6 


A x 


A 2 


As 


10 6 


3.0 


-1.6 


0.9 


1.1 


-0.8 


-2.1 


3.6 


0.5 


-0.5 


0.9 


10 7 


0.95 


-0.1 


0.0 


0.2 


-0.4 


-0.1 


0.0 


0.1 


-0.1 


0.6 



TABLE V: Integrated asymmetries with 10 6 and 10 generated SM events in units of 10 3 . For 
comparison we show in the second column the expected 3a statistical error also in units of 10 -3 . 

CP-odd asymmetries. 

We next turn our attention to the case of strong phases. We show in Table IVJ the 
asymmetries obtained by generating 10 6 events with Madgraph using the cuts in Eqs. [TUJ, [TTJ 
for two different values of /sin Sf. If any of the numbers obtained with /sin Sf = 5 x 



/sin S f 


3 


Ax 


A 2 


^3 


A 4 


A 5 


A 6 


Ax 


A 2 


A 3 


5 x 10~ 4 GeV" 1 


3.0 


0.3 


-2.1 


1.4 


-0.6 


-0.2 


0.3 


1.3 


2.0 


0.9 


5 x 10- 3 GeV" 1 


3.0 


1.3 


-1.0 


-0.2 


-1.6 


2.2 


0.1 


1.2 


-1.3 


1.3 



TABLE VI: Integrated asymmetries with 10 6 generated events in units of 10~ 3 for two different 
values of fsinSf. For comparison we show in the second column the expected 3a statistical error 
also in units of 10 -3 . 

10~ 4 GeV -1 arose from this coupling and not from statistical fluctuations, we would see its 
value increase by a factor of ten for / sin Sf = 5 x 1CF 3 GeV -1 as the asymmetries are linear 
in the anomalous couplings. Table [VI] thus shows, as expected, that strong scattering phases 
cannot fake the truly CP-odd observables we have defined. 

Finally we consider the case of the asymmetries ^1,2,3 as related to CP violation in the 
decay vertex. We show in Table IVHI the results of simulations with 10 7 events using the 
cuts in Eqs. [TUl [TTJ, for two different values of /sin^j. These results appear to show that 
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/sin 4> f 


3 


M 


A 2 


^3 


5 x 10" 4 GeV- 1 


0.95 


2.5 


-0.4 


-1.7 


5 x 1(T 3 GeV" 1 


0.95 


23 


-5.1 


-10.1 



TABLE VII: Integrated asymmetries with 10 6 generated events in units of 10 -3 for two different 
values of / sin </> f . For comparison we show in the second column the expected 3a statistical error 
also in units of 10 -3 . 



these asymmetries are not zero but rather very small. The expected linear scaling with the 
size of the anomalous coupling is not satisfied perfectly leading us to think that statistical 
fluctuations are not completely under control in this case. As stated in the main text, 
however, there are better asymmetries to measure this coupling so we do not pursue the 
issue further. 

To summarize, we have performed the following checks on our numerical results: 

• By generating event samples of different sizes ranging from 10 5 to 10 7 for different 
values of d, and /, we have verified that the asymmetries scale linearly with the 
respective anomalous coupling and that the statistical error scales as the square root 
of the number of events. 

• Using the form factors corresponding to CP violation due to a Higgs boson we have 



reproduced the asymmetry in Ref. 19 



For the case of d and asymmetry A\ we have verified that we are in substantial agree- 
ment with the much more detailed analysis of Ref. [131 ] . 



For the case of strong interaction phases induced by QCD corrections we have checked 



that we are in rough agreement with the estimates in Ref. 
observables. 



14 



which use different 
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